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The current status of our understanding of dilepton production in ultrarelativis- 
tic heavy-ion collisions is discussed with special emphasis on signals from the (ap- 
proach towards) chirally restored and deconfined phases. In particular, recent results 
of the CERN-SPS low-energy runs are compared to model predictions and interpreted. 
Prospects for RHIC experiments are given. 



I. INTRODUCTION 



The rich physics potential of electromagnetic observables (photons and dileptons) has been clearly 
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demonstrated over the past decade of heavy-ion experiments, especially at the CERN-SPS. Although 
the final spectra constitute a superposition of emission from the entire space-time history of a nucleus- 
nucleus collision, the original idea [[lj of identifying signals from (thermalized?!) hot and dense 
phases of strongly interacting matter could be realized. Indeed, an excess of radiation in central 
' Pb(158AGeV)-Au and Pb(158AGeV)-Pb collisions as compared to baseline p-p and p-A measurements 
has consistently been observed, most prominently for electron pairs at low invariant mass [M < 
1 GeV) by CERES/NA45 §, muon pairs at intermediate mass (1.5 GeV < M < 3 GeV) by NA50 § 
and photons at transverse momenta qt > 1.5 GeV by WA98 Q|. ^From the theoretical side it should 
H ' be emphasized that both dilepton and photon emission rates of a thermal medium are based on the 
^ | same quantity, i.e., the (imaginary part of the) correlation function of the electromagnetic (e.m.) 
■ current (or photon selfenergy), Im II cm . Thus a consistent description of both photon and dilepton 
observables within a common approach (coupled with a realistic space-time evolution for a heavy-ion 
collision) is mandatory (see refs. ^|J|] for recent reviews). Whereas real photons are related to the 
M 2 — > + limit of n cm [], the production of timelike photons (^dileptons) requires M 2 > 0. The latter 
^vq , therefore carry additional dynamical information as, e.g., encoded in the low-lying (nonperturbative) 
vector-meson (J p = 1~) resonance excitations (p, uj and </>) of the QCD vacuum. This renders the 
low-mass region an ideal regime to study pertinent in-medium effects in connection with (the approach 
towards) chiral symmetry restoration, as will be further detailed below. At invariant masses beyond 
1.5 GeV, the e.m. correlator becomes perturbative in nature, characterized by a rather structureless 
Oh, continuum with a strength governed by free qq with controllable corrections in the strong coupling a s , 
temperature and density. With an approximately known emissivity the virtue of the probe resides in 
the magnitude of the signal as an indicator of the temperatures reached in the early stages. The same 
idea is the main motivation |?J behind direct photon measurements. There is, however, a subtle but 
important difference: the leading contribution to photon production is of order 0(aa s ) (a = 1/137: 
e.m. coupling constant), whereas the one for (intermediate- mass) dileptons is 0(a 2 ) and therefore 
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under better theoretical control (albeit experimentally suppressed). 



II. ELECTROMAGNETIC RADIATION: HADRON GAS VS. QGP 



A. Production Rate and E.M. Correlator in Free Space 



To leading order in a, the emission rate of lepton pairs l + l (I = e, fi) from thermally equilibrated 
matter per unit 4-volume and 4-momentum reads || 

^jDtherm 2 

~^q~ = f Bose (qo;T) Imn em (M,g;/i B ;T) (1) 



1 The e.m. correlator is also directly relevant to studies of charge fluctuations. The fluctuation content of a 
locally thermalized system is directly proportional to the static spacelike limit of n em . 
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(M 2 = <7q — q 2 ). The theoretical objective is to calculate the (thermal expectation value of the) 
correlation function of two e.m. currents, 

n&(o>>9) = / d 'v * iq ' v m m (y)jonM\o) , (2) 

which, in principle, contains all orders in the strong interaction. In practice, one has to invoke 
approximations starting from either quark or hadronic degrees of freedom. Accordingly, the current 
operator takes the form 



Jem 



J_ 9 'M 4- v2o 

V2 J P 6 J 



3-^ 



hadronic basis 



is7 M s quark basis 



(3) 



ImII om (s) 



In free space, the 'transition' can be inferred from the famous e + e~ — > hadrons cross section: at 
low invariant mass, \fs = M < 1.2 GeV, it is saturated by p, to and <f> mesons (Vector Dominance 
Model = VDM), whereas beyond M ~ 1.5 GeV, the cross section is determined by weakly correlated 
qq pairs with almost no impact from subsequent hadronization (until additional flavor thresholds are 
reached where nonperturbative bound state formation occurs). Thus, the e.m. correlation function 
may be decomposed as 

E (jif) ImDy(s) , s < s dual 

(4) 

( i + ^£) + ...) Nc E (e?)2 , s > SduaJ . 

For the investigation of medium effects through dilepton production, an important question is 
how the spontaneous breaking of chiral symmetry (SBCS) manifests itself in the e.m. (or vector) 
correlator. A direct connection can be established in the isovector (/ = 1) channel: at low mass, 
the pertinent hadronic current is entirely governed by the /c»(770)-meson which has a well-defined 
'chiral partner' in form of the ai(1260)-meson, which itself saturates the low-mass part of the (/ = 1) 
axialvector current. Recent data from r-decays by the ALEPH collaboration |J nicely exhibit this 
feature, cf. Fig. [l]. 
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FIG. 1. ALEPH data for axial-/vector correlators as extracted from r decays 
calculations using p and ai resonances plus perturbative continua. 



The curves are model 



For s < Sdual, the isovector hadronic (axial-) vector correlator can therefore be written as 



ImIPv=j( S ) = 



,IJ P = ir 



(5) 



Im£> 0l (s) - sffr5(s - ml) , IJ P = IV 
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The difference between the two is the manifestation of SBCS in the QCD vacuum. Note the appear- 
ance of the pion pole contribution (the Goldstone boson of the symmetry breaking) with a residue 
given by the pion decay constant. At higher masses both correlators merge into a continuum value 
determined by perturbation theory: 

imn^i = (1 + + . . .) Nc \ , (6) 

i.e., beyond s ~ Sduai effects due to SBCS no longer play a role. 

Chiral symmetry breaking can be further quantified in terms of the Weinberg sum rules ]ic| , which 
are energy- weighted moments of the difference between vector and axialvector correlators, e.g., 

fl = ~ /— (Imn p -ImIT ai ) . (7) 

J 

Kapusta and Shuryak jllj have shown that these also hold in hot and dense matter which provides 
useful constraints on hadronic models. Chiral restoration obviously requires the two correlators to 
degenerate over the entire mass range, i.e., ImiTp(s) = ImiT Ql (s). How this is realized constitutes 
one of the main questions in the context of low-mass dilepton observables and related theoretical 
analyses. In this respect it is fortunate that the isovector channel dominates the strength in the e.m. 
correlator; for two flavors, e.g., the coefficient i in eq. (|^) amounts to ^90% of the total emissivity 
(coefficient | in eq. (^)). 

B. Medium Effects 

Following the decomposition suggested in the previous section, a schematic overview of in-medium 
effects is given in Tab. 0. 

As emphasized before, the low-mass radiation from hadronic matter is well suited to study in- 
medium modifications of light vector-mesons. Their vacuum properties are intimately related to 
the spontaneous chiral symmetry breaking; e.g., My ~ 2Af 9 , with the constituent quark mass 
M q = 400 — 500 MeV driven by the (qq)-condensate. Thus, even below T c , medium modifications 
of vector mesons ought to be considered as precursor phenomena of chiral restoration. Indeed, for 
low-mass lepton-pair emission the overall thermal Bose factor in eq. (Q) induces a relatively moderate 
dependence on temperature, which is largely compensated by the volume increase of the expanding 
and cooling fireball in a heavy-ion collision. Thus, one is rather sensitive to medium effects in the 
spectral distributions themselves. These have been extensively studied in the literature. Starting 
point are typically effective, chirally invariant lagrangians for hadronic interaction vertices with cou- 
pling constants and formfactors determined by experiment. Based on nuclear and finite-temperature 
many-body techniques one evaluates a vector-meson selfenergy £ p = X p7r ,r + EpM + S p s (simi- 
lar for ui and </>) arising from interactions with the surrounding hot (M = tt, K, p, . . . ) and dense 
(B = N, A, A, .. .) hadron gas (HG); £ p7r7r incorporates medium modifications of the free p — > tttt 
decay. The generic outcome of such calculations is a strong increase of |Im£ p | with temperature and 
density, which broadens the p-meson spectral function, Im D p , beyond recognition of any resonance 
structure, cf. Fig. |2|. At comparable densities, nuclear effects prevail over the ones induced by thermal 
pions (roughly speaking, interactions with pions are 'Goldstone'-protected). 



TABLE I. Summary of medium effects in thermal dilepton production. 
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FIG. 2. Real and imaginary part of the p-meson propagator |2fj] at finite temperatures and baryon densities. 

In the invariant-mass window between 1 and 1.5 GeV, the main source of hadronic dilepton pro- 
duction is due to 7rai annihilation This, in fact, can be understood as a pion-induced mixing 
between the axialvector and vector correlator (and vice versa) to lowest order in temperature, char- 
acterized by the parameter e = T 2 /6/ 2 (chiral limit) Jig ]. When extrapolated to high temperatures, 
this 'chiral mixing' leads to a mutual degeneracy of V- and ^-correlators not too far from the critical 
temperature as extracted from lattice gauge theory, T c — 170-190 MeV. 

Another line of treating medium effects employs chiral lagrangians in mean-field approximation 
coupled with density-dependent corrections to masses and coupling constants. E.g., within the so- 
called Brown- Rho scaling scenario jl7]], all hadron masses (except for Goldstone bosons) 'drop' with 
increasing temperature and density. 

Turning to the QGP phase, it turns out that, at low mass, the correct leading-order in a s result 
requires a resummation of thermal propagators accompanied by vertex corrections, which is achieved 
within the so-called Hard Thermal Loop (HTL) framework |l^]. Although the strong coupling con- 
stant g s is not really small for practical applications in the QGP phase, T ~ (1-3) T c , the main quali- 
tative result of a substantial enhancement of the rate towards small M (due to 'Bremsstrahlung'-type 
processes) should be rather robust. At higher masses, M 3> T, naive perturbation theory becomes 
applicable again, with moderate corrections |I§| , |l9| ]. 

Fig. ^ summarizes the various features graphically in terms of the 3-momentum integrated produc- 
tion rate 

5§t " -SCP /§ '"'"^ ™$?Mr.»,.T) (8) 

in the isovector channel, which amounts to a comparison of 7T7t annihilation (or p decays) in the HG 
versus qq annihilation in the QGP phase. Starting from the free 7T7t and qq rates, which are obviously 
very different from each other, the respective in-medium corrections lead to (i) a characteristic low- 
mass enhancement in both the HG and QGP, as well as (ii) p-resonance melting (around M = m p ) 
and chiral mixing (1 GeV< M < 1.5 GeV) in the HG. As a result, the overall emissivities from the 
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hadronic and Quark-Gluon phase look surprisingly similar around T c , which has been interpreted as 
an in-medium reduction the quark-hadron duality scale Sduai Po| • 




FIG. 3. Equilibrium dilepton production rates from hadronic and QGP matter. 



III. DILEPTONS AT CERN-SPS ENERGIES 



A. Comparison to Data 

Under the assumption that the highly excited matter formed in central heavy-ion reactions reaches 
thermal equilibrium at some typical formation time To (~ 1 fm/c at SPS), hydrodynamic or thermal 
fireball simulations are the appropriate framework for a consistent application of the thermal rates 
discussed in the previous section. Based on a fixed entropy per baryon as inferred from produced 
particle abundances, an isentropic thermodynamic trajectory T(/is) consistent with standard hadro- 
chemical freezeout analyses f|l],^2| can be constructed including a phase transition from the QGP 
to HG. An important ingredient in the fireball expansion from chemical freezeout, (T c h, fig) to the 
thermal one, (T t ^, is the conservation of the observed particle multiplicities. This necessitates the 
build-up of finite meson-chemical potentials towards T t h, e.g., /i^" ~ 60-80 MeV at Ei a t, = 158 AGeV. 
The contribution of thermal radiation to observed spectra then takes the form 

dN thermal T f ,,/"., M dR* hermal 



dM 



/C A f rl T?t tier mat 

drV FB (r) d 3 q- dKl+ J- Ace, (9) 



with Vfb the fireball volume, and the factor Acc accounts for experimental acceptance cuts. In 
addition, contributions from long-lived hadron decays after freezeout have to included (the so-called 
'cocktail' 0). 
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FIG. 4. Thermal fireball calculations [|20| for (semi- 
for 40 AGeV compared to CERES/NA45 data " 
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central Pb+Au at 158 AGeV and pertinent predictions 



An example of a thermal fireball calculation plus cocktail is shown in the left panel of Fig. [| 
and compared to CERES/NA45 data from Pb(158 AGeV)+Au collisions. The enhancement over 
the cocktail cannot be described by adding free 7T7t annihilation within the fireball. With the strong 
medium effects as displayed in Fig. || (p-'melting') a reasonable description is obtained. However, 
also the assumption of a dropping p-mass reproduces the data. The QGP contribution is small and 
insensitive to initial temperature and details of the phase transition construction. In the right panel of 
Fig. ^ the predictions of the same approach |p(if are found to be consistent with the measurements at 
lower SPS energies (40 AGeV) as well. In fact, despite the lower pion multiplicities, the calculations 
imply a slightly larger signal in the M ~ 0.4 GeV region due to larger medium effects induced by 
higher baryon densities in the hadronic phase. This is in line with the trend of the data. 
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FIG. 5. Transport calculations |25| for Pb(40 AGeV)+Au; the 158 AGeV data 
orientation only. 



have been included for 



Another example of a low-mass dilepton prediction for the CERES 40 AGeV run is shown in Fig. [| 
in terms of a UrQMD transport calculation p5| |. Closer inspection reveals that the Dalitz decay 
contributions uj — > 7r°e + e _ and 7] — > 7e + e~, which are prevalent around M ~ 0.4 GeV, are a factor 
of 3-4 above the standard hadrochemical cocktail (24] . 

Finally, we display in Fig. [j| a thermal fireball calculation 26 for the intermediate-mass dimuon 
spectra measured by NA50 in central Pb(158 AGeV)+Pb. The 'dual' dilepton production rate based 
on eq. (0) has been folded over the same space-time evolution (modulo centrality) underlying the 
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results of Fig. ffl. As anticipated, the yield is much more sensitive the early phases with a significant 
part of the yield emerging from a QGP with preference for initial temperatures around To ~ 220 MeV. 
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FIG. 6. Fhermal fireball calculations |£6[ compared to NA50 data ||. Upper and lower total yields corre- 
spond to T = 195 MeV and 220 MeV, respectively. 



B. Consistency of Models 

Let us first address the low-mass enhancement observed by CERES/NA45. Under the premise 
that the relative importance of baryon density effects around midrapidity increases towards lower 
collision energies (larger baryon stopping and lower initial temperatures), the indication for a large 
enhancement in the 40 AGeV data suggests that its origin lies in baryon-driven in-mcdium effects 
as predicted within the many-body approach of refs. ]2J|,^0| (unless dramatic nonperturbative effects 
occur in the QGP phase close to T c , which seems unlikely and finds no support in lattice results pTfl). 
At the same time, transport calculations predicted very little effect from higher baryon density p5fl 
and, in their present form, cannot describe the 158 AGeV and 40 AGeV data simultaneously. 

Secondly, the NA50 dimuon enhancement at intermediate mass can be explained by thermal ra- 
diation within the same (fireball) evolution scenario p(| as the CERES data without additional 
parameters, albeit a with higher sensitivity to early (QGP) phases, as discussed above (see also 
refs. P^,Ej|). In particular, an 'anomalous' open-charm enhancement (over the expectation from N- 
N collision scaling) is not required. In fact, also the J/^S data of NA50 J3(J follow within a common 
framework which combines suppression in a QGP and thermal production at T c |32] . 

Finally a remark on thermal photons is in order ||.The WA98 data || seem to require higher initial 
temperatures than the NA50 dimuons. However, photon production (both the hard component from 
N-N collisions as well as the QGP rates) has a nontrivial leading a s -dependence which implies 
relatively larger uncertainties than in the dilepton sector. This equally applies to effects of (nuclear) 
/ct-broadening (Cronin effect), which might be larger than assumed in current calculations and thus 
responsible for (part of) the discrepancy. 



2 Note that the experimental fact that hadron production in e + e~ annihilation follows statistical (thermal) 
model predictions justifies the use of the perturbative ('dual') rate in the intermediate mass region even 
in the hadronic phases of a heavy-ion collision: the thermal heat bath is the same state as produced in 
e + e~ — > hadrons, so that one can use time-reversal invariance to obtain the e + e~ production rate. 
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C. Chiral Restoration? 



Let us briefly elaborate on some recent developments. 

In ref. |}3| the so-called 'vector manifestation' of chiral symmetry restoration has been suggested: 
within the Hidden Local Symmetry framework, the chiral partner of the pion is identified with the 
longitudinal component of the p meson. When applied within a one-loop renormalization group 
evolution towards the symmetry restoration point, the p mass along with and g p go to zero. In 
particular, VDM does not hold at finite T. 

Rather different results are obtained in ref. |34| , where a linear anpai model, consistently con- 
structed to 1-loop order with realistic vacuum properties, has been applied at finite temperature. It 
is found that, even if the vacuum p mass at tree level is entirely given through the (qq) condensate, 
the vanishing of the latter at T c does not si gnifi cantly impact m p (T c ) due to different (leading) tem- 
perature dependencies of the two quantities |16| . The main effect was rather a substantial broadening 
of both the p and a\ spectral functions. This is not inconsistent with the chiral restoration scenario 
put forward in ref. [po[ , where the degeneracy of vector and axailvector correlators is realized through 
the in-medium reduction of the quark-hadron duality threshold. 

The relative importance of baryonic effects in 'melting' the p resonance calls for a better un- 
derstanding of their relation to chiral symmetry. A possible scheme is illustrated by the diagrams 
in Fig. |?], which combines pionic S'-wave excitations connecting chiral partners with pionic P-wave 
excitations representing hadronic resonance physics. 
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FIG. 7. Interaction scheme combining chiral and resonance excitations. 



In particular, the analogy between the mesonic (left panel) and baryonic (right panel) sector iden- 
tifies the prominent role of the pN S-w&ve resonances /V(1520) and A(1700). 



IV. PREDICTIONS FOR RHIC 



The beginning of the collider era through the first operation of BNL-RHIC has opened a new 
energy frontier in heavy- ion physics. First dilepton data are expected from the second run of Au+Au 
(completed in Nov. 2001) at the full energy of ^fs = 200 AGeV. 

Predictions for e + e~ spectra in central Au+Au are summarized in Fig. || for the invariant-mass 
range from to 2.5 GeV. The thermal contribution has been evaluated within a fireball model 
with realistic values for the charged particle multiplicity and an estimated p/p-ratio of ~75%. 
The radiation from the QGP phase (with initial temperature To ~ 380 MeV) dominates the thermal 
yield at masses M > 1.5 GeV. Its detectability will critically depend on the contribution from 
correlated open-charm decays, which in Fig. |§| has been taken from a PHENIX event generator |j37fl 
which is based on an extrapolation of N-N collisions (using PYTHIA) at lower energies. If the 
spectrum of the c-quarks experiences significant softening (e.g., through energy loss in the QGP), 
the associated dilepton spectrum may be severely suppressed above M = 2 GeV or so pM), opening 
the window for the QGP signal. In the low-mass region thermal radiation mostly originates from 
the hadronic phase and compares favorably to both open-charm and the hadronic decay cocktail 
after freezeout |37j (here, subtraction of the combinatorial background will be the main experimental 
problem). In-medium modifications again 'melt' the p-resonance which drives the enhancement below 
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M ~ 0.6 GeV 0. Similar effects lead to an in-medium signal from w-decays with an average width of 
~ 50 MeV, almost seven times its vacuum value (less pronounced for the 4>). 
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FIG. 8. Dilepton spectra from central Au+Au collisions at full RHIC energy. 



V. CONCLUSIONS 

Driven by exciting data, the theoretical and phenomcnological analyses of dilepton production in 
high-energy reactions of heavy nuclei have undergone continuous progress with increasing consensus 
among most approaches. As for the low-mass region, the CERES/NA45 data at full SPS energy 
(158 AGeV) require strong medium modifications of the p-meson indicating that one is indeed probing 
strong interaction matter in the vicinity of T c (with small contributions from the QGP phase itself). 
The new 40 AGeV data support the importance of baryon-driven effects, while their relation to chiral 
symmetry (especially for resonances) remains to be better understood, 'p-melting' with 'quark- 
hadron duality' towards T c remains a viable scenario of chiral restoration. At intermediate masses, 
the NA50 enhancement consistently emerges from the same thermal source, but with significantly 
larger sensitivity to the QGP phase pointing at initial temperatures To > 200 MeV. New exciting 
insights can be expected soon from PHENIX at RHIC, HADES at SIS and farther into the future 
from ALICE at LHC. This raises the hope of a systematic understanding of dilepton emissivities 
across the QCD phase diagram. 
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